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The adsorption or adhesion of large particles (proteins, colloids, cells, . . . ) at 
the liquid-solid interface plays an important role in many diverse applications. 



'^ ■ Despite the apparent complexity of the process, two features are particularly 

^ . important: 1) the adsorption is often irreversible on experimental time scales 

O I and 2) the adsorption rate is limited by geometric blockage from previously 
adsorbed particles. A coarse-grained description that encompasses these two 

. properties is provided by sequential adsorption models whose simplest exam- 

QQ ■ pie is the random sequential adsorption (RSA) process. In this article, we 

Cn , review the theoretical formalism and tools that allow the systematic study 

\^ [ of kinetic and structural aspects of these sequential adsorption models. We 

^D ' also show how the reference RSA model may be generalized to account for a 



variety of experimental features including particle anisotropy, polydispersity, 
bulk diffusive transport, gravitational effects, surface-induced conformational 
and orientational change, desorption, and multilayer formation. In all cases, 
the significant theoretical results are presented and their accuracy (compared 
'^ I to computer simulation) and applicability (compared to experiment) are dis- 

Q ■ cussed. 

O 

> '■ 
•rH , I. INTRODUCTION 

X 

c^ I The adsorption of large particles (proteins, viruses, bacteria, colloids, macromolecules) 

at the liquid-solid interface has received considerable theoretical and experimental attention 
in recent years. Among the many issues pertaining to this problem, we address in this 
review the consequences on the adsorption kinetics and the adsorbed layer structure of two 
key features: i) the absence of reversibility and ii) the deceleration of adsorption due to 
surface exclusion from previously adsorbed particles fTj-Q. We choose as a starting point a 
mesoscopic approach, where atomic-level detail is coarse-grained and incorporated as a set of 
effective parameters such as rate constants, particle shapes and sizes, effective interparticle 
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interaction parameters, etc. . . The two essential features listed above are well accounted for 
in sequential adsorption models [0,§; these will be the central topic of this review. 

The outline of our review is as follows. In section |I|, we introduce the simplest of 
sequential adsorption processes, the random sequential adsorption model (RSA), and its 
one-dimensional version, the car parking problem. Since these processes result in structures 
that are generally far from equilibrium, the well-developed methods of liquid-state statistical 
mechanics are not directly applicable. We present, in section |ITT|, the formalism and prin- 
cipal theoretical tools that we have developed for describing the kinetics of the adsorption 
processes and the structure of the adsorbed layers. The next sections are devoted to exten- 
sions of these methods to processes involving additional physical features. In section 0, 
we consider the adsorption of nonspherical particles and of mixtures and in section 0, we 
address bulk transport issues such as particle diffusion and the influence of the gravitational 
field. In section ^I[ we treat surface-induced events, like conformational and orientational 
changes of the particles and in section |V11| , we introduce partial reversibility produced by 
desorption. Finally, we discuss multilayer formation in section [VIII| . 

II. RSA: A REFERENCE MODEL 

The random sequential adsorption model, the prototype for sequential addition processes, 
is a stochastic process in which "hard" particles are added sequentially to a D-dimensional 
volume at random positions with the condition that no trial particle can overlap previously 
inserted ones. 



A. The car parking problem 

The one-dimensional version of the model, known as the car parking problem, was first 
introduced by A. Renyi |^, a Hungarian mathematician. Consider an infinite line, assumed 
empty at t = 0. Hard rods of length a are dropped randomly and sequentially at rate ka onto 
the line, and are adsorbed only if they do not overlap previously adsorbed rods. Otherwise, 
they are rejected. If p{t) denotes the number density of particles on the line at time t, the 
kinetics of this process is governed by the equation 

^ = kMt), (1) 

where $(t), the insertion probability at time t, is the fraction of the substrate that is available 
for the insertion of a new particle. To calculate this quantity, it is convenient to introduce the 
gap distribution function G{h,t), which is defined so that G{h,t)dh represents the density 
of voids of length between h and h + dh at time t. For a given void of length h, the available 
length for inserting a new particle is h — a, and therefore the available line function $(t) is 
merely the sum over the number density of available intervals, i.e. G{h,t): 

<l>(t) = / dh{h-a)Gih,t). (2) 

Since each interval corresponds to one particle, the number density of particles p{t) can be 
expressed as 



p{t) = / dhG{h,t), (3) 

Jo 

whereas the uncovered hne is related to G{h, t) by 

POD 

l-p{t)a= / dhhG{h,t). (4) 



These two equations, Eqs. (^ |), represents two sum rules for the gap distribution function. 
During the process, the gap distribution function G{h, t) evolves as 



dG{h,t) 
d{kat) 



-H{h-a){h-a)G{h,t) + 2 dh!G{h',t), (5) 

Jh+a- 



where H{x) is the unit step function. The first term of the right-hand side of Eq. @ 
(destruction term) corresponds to the insertion of a particle within the gap of length h (for 
h > a), whereas the second term (creation term) corresponds to the insertion of particle in 
a gap of length h' > h + a. The factor 2 is due to the two possibilities of creating a length 
h from a larger interval h'. It is worth noting that the time evolution of G{h,t) is entirely 
determined by intervals larger than h. We now have a closed set of equations, which results 
from the fact that the adsorption of a particle in one gap has no effect on other gaps. The 



above equations can be solved by introducing the ansatz [10 



G{h,t) = ^^!^exp{-Uh~a)t), (6) 

which leads to 

F (t) = t' exp ( - 2 r du ^—^\ . (7) 

Inserting Eqs. (P) and (|^) , one obtains G{h, t) ioi h > a and integrating Eq. (|^) with the 
solution of G{h, t) for h > a finally gives G{h, t) for < h < a, 

2 /■'="°'* F(v) 

G{h,t) = — duexp{-uh/a)^-^. (8) 
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The three equations, (|I]), (^ or (^, all lead to the same result for the number density p(t), 

pit) = -l d„exp(-2y^d.^— j. (9) 

which was first derived by Renyi 0. 

A first nontrivial property of the model is that the process reaches a "jamming limit" 
(when t -^ oo) at which the density saturates at a value Poo^" = 0.7476. . .; this value is 
significantly lower than the closed-packed density (pooC" = 1) that is expected when the 
surface is allowed to restructure between adsorption steps. Moreover, it is easy to show that 
the jamming limit depends upon the initial configuration, here an empty line. In contrast, 
the final state of an equilibrium system is determined solely by the chemical potential and 
has no memory of the initial state. The long-time kinetics, during which a small number of 



adsorption events occurs, can be obtained from Eq. (j^) and are seen to display a power-law 
behavior. 

where 7 is the Euler constant. 

The structure of configurations generated by this irreversible process has several notice- 
able properties. At saturation, the gap distribution function shows a logarithmic divergence 
at contact, h ^ 0, 

G(/i,oo)~ -e-27in(/i/a). (11) 

The correlations between pairs of particles are extremely weak at long distances. The pair 
correlation function behaves at long distances as 

1 f 2 V^" 

g{r) - 1 (X — -- , (12) 

1 [r/cr) \[n r/a J 

where r(a;) is the Gamma function, i.e., it decays super-exponentially in contrast with the 
characteristic exponential decay of equilibrium systems |Tl| . 



B. Two-dimensional sequential adsorption 

The RSA model can be studied in two dimensions, which represents the physical di- 
mension for adsorption problems. Since the RSA process generates random disordered con- 
figurations of "hard particles" , a statistical geometric approach is useful for describing the 
system. For a given time t, or equivalently a given density p, it is convenient to introduce 
the set of ra— particle density functions {p("^)(ri, r2, ■ ■ ■ , r„, t)}, where ri,r2,---,r„ denote 
the positions of the n particles, that characterize completely the configuration of particles. 

The time evolution of the density is still given by Eq. (|I|), where $ denotes the available 
fraction of the surface for the insertion of a new particle. When the particles are hard disks, 
$ has a simple geometrical interpretation: It represents the probability of finding a circular 
cavity of diameter equal to twice the particle diameter, that is centered at ri, and is free of 
particle centers ||T2|-[l^ 



°^ 1 /■ f 
•^(rl; p) = ^-, ■■■ dv2--- rfr(,+i)/i2 ■ ■ ■ /i(s+i)p(')(r2, ■ ■ ■ , r(,+i); p), (13) 

s=0 ^- •' •' 

where * denotes a cavity (as defined above) and fij = /(|rj — r^l) is the Mayer function for 
non overlapping hard particles: 

^• = 0, \vi-rj\ >a, (14) 

fij = -1> ki-rj| < a, (15) 

where a is the diameter of the hard spherical particles. (Note that because of the macroscopic 
uniformity of the system, <l>(r^;p) is independent of the position ri, i.e., <l'(r^;p) = '^(p)-) 



The terms of the sum in the right-hand side of Eq. (0) can be interpreted geometrically : the 
first term, which is equal to one, corresponds to the certain acceptance of particles added to 
an empty surface. The second (negative) term subtracts the fraction of the surface occupied 
by exclusion disks in which no new particle center can be placed. The third (positive) term 
corrects for the fact that two exclusion disks may overlap, etc... For hard-core interactions, 
the sum is always finite since the number of non-overlapping particles located at a distance 
less than the particle diameter from a given point is finite. For instance, for hard disks in 
two dimensions, only six first terms have a nonzero contribution to Eq. (pIS]). 

Higher order density functions can be systematically defined. For instance, $'^^)(r^, r2; p) 
represents the density of finding one (unspecified) particle at the point Y2 and a cavity of 
diameter 2cr free of particle centers at the point ri. More generally, one can introduce 
$*^"'^(r^, r2, ...Tn',p) as the probability density of finding n — 1 (unspecified) particles at the 
points r2, r3, ..., r„ and a cavity of diameter 2a free of particle centers at the point ri. 



$(")(r*,r2,---r„;p) = n(l + /: 



ii> 



i=2 



ij2^ j ■■■ j drn+i ■ ■ ■ dv(^ri.+s)fi(n+i) ■ ■ ■ fi(n+s)P^''^' ^^ (rg, ■ ■ ■ , r(„+,); p) j . (16) 

For a RSA process, the insertion of a new particle amounts to finding a cavity of diameter 
2cr for a given set of previously adsorbed particles. The kinetic equation for the density, 
Eq. (|l]), can be generalized to higher order density functions, 

<9pn(ri,r2, .•.,r„;p(t)) " . s 

— — = 2.$^ ^(ri,.,ri_i,r^,ri+i,..r„;p), (17) 



i=l 



or, by using Eq. (|l|), 

c}p„(ri,r2,...,r„;p) 1 



X:$^"nri,.,r._i,r*,r,+i,..r„;p). (18) 



dp $(ri;p)^ 

The exact and infinite hierarchy of equations cannot be solved (nor is it possible in the 
equilibrium case), but approximate solutions and computer simulations (both to be detailed 
below) provide a good description of the process. 

The main features observed in the one- dimensional car parking problem persist in higher 
dimensions: existence of a jamming limit at which the density of adsorbed particles saturates 



(in 2D, ^00 = Poo^ — 0.547 [jl5|,|T6|), slow kinetics when approaching jamming (in 2D, 9oo — 
9{t) ~ t~^/^ EH), logarithmic divergence at contact of the pair correlations at saturation, 
configurations of adsorbed particles that differ from those characteristic of an equilibrated 
system. The qualitative and quantitative relevance of two-dimensional RSA for describing 
the adsorption of macromolecules at liquid-solid interfaces has been demonstrated in several 



experiments on proteins [[l|,p!7| and colloidal particles P,|6|,p!8|-p0| . 



III. THEORETICAL FORMALISM AND TOOLS 



A. Statistical mechanics of sequential adsorption systems 



Consider a system of A^ particles in a D- dimensional volume V (in practice, adsorption in- 
volves D = 2). For simplicity, we assume that the particles are identical and interact through 
a pairwise-additive spherically-symmetric potential. Generalization to mixtures, nonspher- 
ical particles and non pairwise-additive potentials is, conceptually at least, straightforward 
(see below). If the system of particles is at equilibrium at a given temperature T, the proba- 
bility density, p[^\{ri, r2, ■ • -rTv}) associated with finding the A^ particles in a configuration 
with one (unspecified) particle whose center is an infinitesimal volume element around posi- 
tion ri,. . . ,one (unspecified) particle whose center is an infinitesimal volume element around 
position tn is given by the Gibbs distribution for the canonical ensemble. 



Pir({ri,r2,---M) 



ZNiV,(3) 
l/ksT, ks is the Boltzmann constant, f (|rj — 



where /3 

Zn{V,i3) is the configurational integral defined by 



(19) 
is the pair potential, and 



Zjv(V, /3) = / • • • / dridr2 ■ ■ ■ dT^e 



-PH. 



<i<j<N 



v{\ri-rj\) 



(20) 



All the powerful tools of equilibrium statistical mechanics and the connection to thermo- 
dynamics derive from the Gibbs distribution. Consider now a situation in which the con- 
figuration of particles does not correspond to thermal equilibrium, but is generated by a 
sequential adsorption process. More specifically, we consider a "cooperative sequential ad- 
sorption" (CSA) that is a simple generalization of the random sequential adsorption process 
introduced above. In this process, the rate of adsorption of a new particle in a given con- 
figuration of pre-adsorbed particles is proportional to a Boltzmann factor involving the 
interaction energy between the new particle at the chosen position and all preadsorbed par- 
ticles. If {ri, r2, ■ ■ ■ r„} denotes the position of the n preadsorbed particles, adsorption of a 
particle at position r„+i is given by 



/c(r„+i|{ri,r2,---r„}) = kaexp 



-pY^^d^n+i-r-j 



(21) 



where ka is the rate of adsorption of a particle in an empty D-dimensional volume. What 
is the probability of finding a given configuration with a total number of N particles at 
{ri,r2, ■ ■ -riy} which is generated by this cooperative sequential adsorption, irrespective of 
the time required? It is clear that the intrinsic irreversibility, resulting from the fact that 
once a particle is adsorbed it can no longer move or desorb, creates a strong memory effect. 
As a result, one must take into account the order in which the positions ri,r2, ■ ■ -tn have 
been filled by particle centers and consider the A^! sequences of insertion that lead to a 
given configuration of A^ particles at {ri, r2, ■ ■ ■ rjv}. With the rate given by Eq. (^) and 
the definition e(r„+i|{ri, r2, ■ ■ -r^}) = A;(rn+i|{ri, r2, ■ ■ ■ r„})/fca, one obtains the following 
probability density 

pS.({r„r,,...r„)} = 



-"E. 



<i<j<N 



v{\ri 



N\ 



— En 

iTS5jv 4 = 1 



(ir^(i)e(r^(i)|r^(i). 



[■(T{i-1) 



}) 



-1 



(22) 



where the sum o G S'at is over all permutations of {1,2, ■■■,A^} and where we have 
used that Ilili e({r<^(i)|r^(i), r^(2), • ■ -r^Ci-i)}) = exp (-/?X;i<i<j<n^(ki - rj|)). (One eas- 
ily checks that the normalization condition J ■ ■ ■ J dridr2 ■ ■ ■ dr]yp}jgj^{{ri, r2, ■ ■ ■ fAf)}) = 1 
is satisfied). Eq. (^) is the CSA counterpart of the Gibbs distribution for equilibrium 
systems, Eq. (|19|). In a thermally equilibrated system, Eq. (|I9|) implies that all config- 
urations having the same energy are equiprobable. On the other hand, in a sequential 
adsorption process there is a bias that favors some configurations over others among all 
those that have the same energy: this results from the fact that in Eq. (|22D , the factor 

EaGSArll^i /c?r^(i)e(ra(i)|r<T(i),---ra(i--i)}) depends on the positions ri,r2,---rjv (for 



A^ > 2). This feature has been already noted by Widom [^] in the example of three 
hard spheres added to a large volume. The seemingly formal difference between the equi- 
librium Gibbs distribution, Eq. (p!9D, and the CSA distribution, Eq. (p2D, has important 
consequences. Some of the dramatic ones are that the CSA systems are characterized by the 
existence of a jamming limit where the density of adsorbed particles (for particles possessing 
a hard core) reaches a saturation value that is significantly less than the optimum filling 
(as in RSA, cf section |1[) and the absence of thermodynamic phase transitions^. This latter 
feature results from the infinite memory of the adsorption process and is illustrated in Fig |l|. 
In the figure, we compare two typical configurations of hard nonspherical two dimensional 
particles (disco-rectangles , i.e. a planar sphero cylinder, with an aspect ratio of 15) on a 
plane at the same surface density (or coverage), one being an equilibrium configuration, the 
other being generated by RSA. At the chosen coverage of 0.445, the equilibrium configu- 
ration is in a nematic phase with long range orientational order (the transition from the 
disordered to nematic occurs for 6 ~ 0.4), whereas the RSA configuration is at saturation 
in a disordered state. 

The existence of a well-defined probability density for configurations of N particles pro- 
duced by sequential adsorption in a D-dimensional volume \^ at a given temperature T, as 
in Eq. (P2|), makes a "statistical mechanical" approach possible, but the complicated form 
of the probability density imposes a nontrivial generalization of the equilibrium formalism. 

The kinetics of a sequential adsorption process, as well as the characteristics of the 
configurations of adsorbed particles, can be obtained from a knowledge of the grand ensemble 
analog of the "canonical" probability density introduced above, p*^^^({ri, r2, • • -rAr}; t), which 
is associated with finding a configuration of A^ adsorbed particles positioned at {ri, r2, ■ ■ ■ r^} 
at time t. However, the expression of this quantity is quite intricate, and it is convenient to 
derive directly kinetic equations describing the time evolution of the infinite set of n-particle 
densities, p'^")({ri, r2, ■ ■ ■ r„)}; t). (p'^")({ri, r2, ■ ■ -Tn}; t) is associated with the probability 
of simultaneously finding an (unspecified) particle at ri, another at r2,. . . , and one at r„, 
irrespective of the remaining particles and is thus obtained from the p^^''({ri, r2, ■ ■ ■ r„}; t)'s 
for A^ > ra by integrating over the positions oi N — n remaining particles. This derivation 
has been done in section [11 B| for the RSA of hard particles by means of statistical geometric 



arguments. The result can be generalized for all sorts of cooperative sequential adsorption 
processes (as well as processes involving desorption and other mechanisms, see the next 



^Geometric phase transitions such a percolation phenomena are still possible, see section VB 



section). Consider a cooperative sequential adsorption similar to that described above with 



an adsorption rate that generalizes Eq. (21), namely 



fc(n+l|{l,2,---n}) = A;„e(n + l|{l,2,---n}) 

= fc,exp[-/?M/(n+l|{l,2,---n})] (23) 

where j is a short-hand notation for center-of-mass position r^ and orientation 6j and 
iy(n + 1|{1, 2, • ■ -n}) is the interaction energy between an incoming particle at position 
Tn+i with orientation 6n+i and n preadsorbed particles characterized by {1, 2, ■ • • n}. There 
is no need here to assume pairwise additivity or spherical symmetry for the interaction 
potential. 

The time evolution of the 1— particle density, which, for macroscopically homogeneous 
(i.e. both uniform and isotropic) systems considered here, is equal to the average density of 
adsorbed particles p, is governed by the following kinetic equation: 

where $(l*;t) is the probability of inserting a new particle at the position ri with an 
orientation 0i. (Note that because of the homogeneity of the system, one can average as 
well over ri and ^i). Introducing the Ursell (or "cluster") functions associated with the 
e(n + l|{l,2,---n})'s, 

/n+l|{l,2,...n} = E (-1)'^-"^ E e(n + l|{ii, i2, ■ ■ ■ i^}), (25) 

"1=0 {ii,i2,...,im}C{l,2,...n} 

where the second sum in the right-hand side of the above equation is over all distinct m- 
tuplets chosen from {1, 2, . . . , n}, and recalling that for sequential adsorption processes there 
is a one-to-one mapping between reduced time, kat, and density, p, one generalizes Eq. (p!3D 
as follows: 

$(1*; p) = Y^- ... d2...d{s + l)/i|{2,..(s+i)}p(^n{2, ■ ■ • (s + 1)}; p), (26) 

where the integral involves integration over both center-of-mass positions and orientations. 
For a pairwise additive potential, it is easy to show that the Ursell function is just a product 
of conventional Mayer functions fij = e~^^^^ — 1, /i|{2,---{n+i)} = 11^=1 fi{j+i)ji cind that 
Eq. (|26| ) reduces to Eq. (|1^) for simple RSA. For higher-order particle densities, generalizing 
Eq. ([ISD leads in a similar way to 

c)pW(1.2....n;,) ^^_" i.ir i+1 n.p) (27) 

with 

$('^)(r , 2, ...n; p) = E -y / • • • / rf(n + 1) . . . rf(n + s)/(\^^-(„"i)...(„+,)„pW(2, ...(n + s); p), 

(28) 



where the /i||(n+i)...(n+s)|'s are generahzed Ursell functions defined by 



s 



/i|«n?i)...(„+s)}= E(-ir""^ E e(l|{2,---n,ii,...,U) (29) 

»"=0 {n,...im}C{n+l,...,n+s} 

Again, for pairwise additive potentials, the generahzed Urseh functions reduce to prod- 
ucts of Mayer functions. (In section |V B| , we shah consider a process with non-pairwise 
additive potentials.) 

The above hierarchy of equations for the n-particle densities is the counterpart for coop- 
erative sequential adsorption processes of the Kirkwood-Salsburg hierarchy for equilibrium 
systems. The latter reads 

^^"^(^'^'•••"'/^) =^_y^M(l i-liM+l n-p) 
n$("ni*,2,...,n;p) 



$(l*;p) 



(30) 



where the $(")'s are given by the same expression as in Eqs. ( PU] ) and (P^). The com- 
parison of Eqs. (|27|) and (|30|) illustrates again the difference between CSA and equilibrium 
systems. 

The Kirkwood-Salsburg hierarchy for CSA processes, Eqs. (p^)- (pHf), is a good starting 
point for generating diagrammatic expansions ||l4|. These can be used to formulate various 



approximate descriptions, such as low-density expansions for the adsorption probabilities and 
integral equations for the pair correlations, that parallel those used in equilibrium statistical 
mechanics of fluids. In the case of pairwise additive potentials, diagrammatic expansions are 
also conveniently derived by using the "replica trick" developed for spin glass models p2| . 



The trick relates the probability distribution for CSA and RSA conflgurations of A^ particles 
to the Gibbs distribution of a flctitious equilibrium system with n2 replicas of particle 2, n^ 
replicas of particle 3,. . . , n^v replicas of particle A^, in the limit n2 — ^ 0,n3 -^ 0,n7v -^ 0. pH 



B. Approximation scheme 

There is, unfortunately, no exact solution to the sequential adsorption problem when 
the dimensionality D of the substrate is two or more, (the one-dimensional case will be 
discussed below). To describe the kinetics of the adsorption process and of the structure 
of the adsorbed conflguration, one can use computer simulation (see below). In addition, 
it is also fruitful to develop approximate treatments that allow one to study a wide range 
of control parameters and a variety of interactions, external flelds, and mechanisms. To do 
so, it is instructive to consider conflgurations of hard disks generated by RSA at different 
stages of the process. This is illustrated in Fig. ^-c. Also shown in white on the flgures is 
the region that is available for the insertion of a new disk center and whose area gives the 
adsorption probability (correspondingly, the region in gray is that part of the surface which is 
excluded to the center of a new disk, over and above the region covered by the adsorbed disks 
themselves that is shown in black, cf section PB| ). At low surface coverage, corresponding 
to short times, the available surface percolates through the whole space and adsorption is 

9 



locally prohibited by groups involving only a small number of preadsorbed disks. When the 
coverage increases, the available surface diminishes and no longer percolates. Finally, close 
to the jamming limit (the saturation coverage is 9oo — 0.547), the available surface reduces 
to small, isolated regions that can accommodate only one additional disk. This suggests 
that different approximations must be used for describing the different stages of the process. 
Consider first the last-stage, long-time, regime that corresponds to the asymptotic ap- 
proach to the jamming limit. Note that the existence of such a jamming limit is a feature 
of all sequential adsorption processes, be they random or cooperative, provided the adsorb- 
ing particles possess a "hard core" that results in an exclusion or blocking effect. When 
the surface coverage is sufficiently high, this hard-core interaction dominates all additional 
"soft" potentials describing for instance dispersion and electrostatic interactions. (The case 
of "singular" additional potentials such as those encountered for systems of sticky disks and 
in the description of generalized ballistic deposition (see section [VB| ) is different). In the 
asymptotic regime near the jamming limit, the surface available for adding new particles is 
composed, as shown above, of small isolated "targets" , whose number exactly provides the 
number of particles that will still be adsorbed before reaching saturation. Due to the irre- 
versibility and the infinite memory of the process, the asymptotic kinetics is simply related 
to the rate of disappearance and the statistics of the targets. The standard arguments for 



the RSA of spherical particles have been put forward by Pomeau [24| and Swendsen [25 
Consider a typical target at time t that is characterized by a small linear scale h so that the 
area available for the center of a new disk goes as h"^ in 2 dimensions (and h^ for general 
dimension D) (see Fig. ^. The rate of disappearance of such a target is proportional to its 
surface area, hence to h?-^ as a consequence, the number density of targets characterized by 
a linear size /i, n(/i, t) decays exponentially in time according to 

n{h, t) ~ n(/i, t,)e-^-^^^''^'-'^\ (31) 

where tc is a crossover time that marks the beginning of the asymptotic regime and A is 
a mean shape factor (besides /i, the geometric characteristics of the targets are irrelevant). 
Making the reasonable assumption that ra(/i, t) goes to a nonzero value n(0, t^) when /;, — > 0, 
one then obtains that the difference in density of adsorbed particles between time t and 
saturation {t — > cx)) is given by 

p(oo) — p(t) = / dhn{h,t) 
Jo 

~ n(0, te) r^ dh e-^-^^""' ~ t-i/^, (32) 

Jo 

where he is an (unimportant) upper cutoff. The approach to saturation is thus very slow 
in RSA processes and is described by a power law in time (when D = 1, one recovers the 
exact t~^ result, cf section [11 A| ). The same line of reasoning permits one to show that the 
pair correlation function at saturation has a logarithmic divergence when two disks (or, 
more generally, two Z)— dimensional spherical particles) approach contact, a result that is 
also obtained exactly in one dimension (cf section |1I A| ). The arguments in terms of small, 
isolated targets can be generalized to all CSA and RSA processes involving particles with 
a hard core. The key property is the rate of disappearance of targets when their linear 
size h goes to zero. If this rate goes to a nonzero limit, as in the case of generalized 
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ballistic deposition (which corresponds to singular sticky potentials, cf section |VB| ) and of a 



mixture involving point-like particles ^Q\, the approach to saturation is fast and is essentially 



exponential. If the rate goes to zero, generally as /i where J\f is the number of degrees of 
freedom for the adsorbing species, the approach is again a power-law, t~^^ : for instance, 
in two dimensions. A/" = 3 for non spherical (unoriented) particles p7|j28|] as well as for a 



continuously polydisperse mixture of disks p9|, whereas A/" = -D = 2 for a monodisperse 
system of disks. Note that the above results are valid for irreversible adsorption in which 
once adsorbed, the particles stay fixed. A quite different behavior is observed when particles 
are allowed to desorb from the surface (see section |V11.B| ). 

In the low and intermediate density regime (Figs. ^ and ^), one expects that the ap- 
proximation employed in the statistical mechanical theory of equilibrium fluids can, mutatis 
mutandis, be useful as well for describing systems of particles produced by sequential adsorp- 
tion. In particular, this is the case for low-density expansions, somewhat analogous to the 
standard virial expansion for the imperfect gas |^|. As mentioned above, such expansions 



can be generated either order by order from the Kirkwood-Salsburg-like hierarchy or directly 
from the diagrammatic series when they are available (e.g., for pairwise additive potentials). 
For instance, the adsorption probability $rsa(^) for the RSA of hard disks onto a plane is 
given by 

where a is the hard-disk diameter and 6 = {tkj'^ /A)p is the surface coverage. As noted before 
(section |IIIA ), this expansion coincides with its equilibrium counterpart at order 6'^, but 



departs from it at the next order (in general, ^rsa < ^eq)- Similar formulas have been 
derived for the adsorption of hard convex objects of various shapes, for mixtures, and for 
sequential processes involving cooperative effects and competing mechanisms (see below). 
For cases involving "soft" (i.e continuous in the configurational variables) interaction poten- 
tials in addition to hard cores, the computation of the coefficients is always numerical since 
it requires more than geometric considerations. However, one can again fruitfully borrow 
ideas that have proven efficient in equilibrium statistical mechanics [^: weak attractive 



potentials can be treated as perturbations and spherically-symmetric soft repulsive poten- 
tials can be approximately replaced by a hard-core potential with a suitably chosen effective 
diameter that depends on both density and temperature P,pT[]. 

Even for hard-core exclusion interactions, calculating the coefficients rapidly becomes 
intractable as the order in density increases. For hard disks, the fourth-order term has 
been obtained and for parallel hard squares, the sixth-order term, but these are exceptions, 
and in practice only coefficients up to p^ (or 6^) have been calculated analytically. This 
unfortunately restricts the range of density over which the truncated expansion is accurate. 
It is possible, however, to enhance the convergence of the low-density expansion by building 
interpolation formulas that incorporate the known asymptotic behavior. For instance, in 
the case of the RSA of hard disks, the result p(oo) — p(t) ~ t^^/^ converts, through the 
use of the kinetic equation dp/dt = <l>(p(t)), to <l>(p) ~ (^(oo) — p)^ when p —^ p{oo). By 
employing standard techniques such as Fade approximants ||3^-B4[, e.g.. 



Z^n=0 '^nP 
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where the {p + 1) coefficients a„ are determined by expanding in coverage (or density) and 
matching term by term to the exact expansion, Eq. (|33|) , one obtains an improved description 
of the kinetics over the whole process. In particular, expressions such as Eq. ( p^ provide 
estimates of the saturation coverage ^(oo). The procedure works very well for the RSA of 
hard D— dimensional spheres and D— dimensional parallel hard cubes: cf Table I. For more 
complicated shapes and more complicated cooperative processes (see below), the results are 
fair, but not as good. This trend can usually be rationalized by invoking the appearance 
of an additional regime in the process (see, for instance, the case of elongated nonspherical 
particles with random orientations discussed in section |1V A]) . 



The approximation scheme just described focuses on the kinetics of the process at the 
level of the density of adsorbed particles (or surface coverage). However, valuable insight on 
the adsorption process is gained by a proper description of the spatial correlations among 
adsorbed particles. Actually, studying the structure of the adsorbed layer can facilitate 
the discrimination between different hypothesized model descriptions of real adsorption sit- 
uations [ll|,|6|,|20[| . In this regard, it is interesting to develop approximate treatments for 



obtaining the pair correlation function in configurations of particles generated by sequential 
adsorption. As in equilibrium statistical mechanics, the most obvious route is to derive 
integral equations for the pair correlations, and this can be done by making approximations 
in the diagrammatic expansions, as obtained for instance from the Kirkwood-Salsburg-like 
hierarchy presented above. As an example, the analog of the celebrated Percus-Yevick in- 
tegral equation for Z)— dimensional equilibrium hard-sphere fluids [^ has been studied for 



the RSA case ^5[ . The result for the pair correlation function for an adsorbed layer of hard 



disks is illustrated in Fig ^. The description is good at intermediate surface coverage, but, 
as expected, it deteriorates in the asymptotic regime near saturation; the approximation 
fails to reproduce the logarithmic divergence at contact (see section [11 Bj), but it captures 



the rapid, super-exponential decay of the spatial correlations with distance that is typical 
of sequential adsorption processes [§|,|n|- 

C. Usefulness of one-dimensional models 

The vast majority of real situations in which large particles (proteins, cells, colloids,. . . ) 
adsorb at a liquid/solid interface correspond to a two-dimensional substrate (that, for sake 
of simplicity, we consider flat and uniform in this whole review) . Why therefore should one 
study one-dimensional models? Equilibrium statistical mechanics seems to teach us that 
one-dimensional systems have very little relevance to physical two- and three-dimensional 
fluids; indeed, no thermodynamic phase transitions are possible in one dimension whereas 
they are a key phenomenon of real fluids. The situation is, however, different for sequential 
adsorption processes. First, as already stressed, there is no thermodynamic phase tran- 
sition in the adsorbed layer, regardless of its dimensionality. Second, the main features 
of the phenomenology of sequential adsorption processes, such as the existence of a jam- 
ming limit and of a well-characterized asymptotic approach toward saturation, the short 
(super-exponentially) decay of the spatial correlations among adsorbed particles, and many 
other properties that may depend upon the specific process under study (see the following 
sections), are independent of space dimensionality. One-dimensional sequential adsorption 
models already have the ingredients, partial or total irreversibility with the associated mem- 
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ory effect, exclusion or blocking phenomena, that make their behavior nontrivial and qual- 
itatively similar to that observed in higher dimensions. There are of course quantitative 
changes. For instance, a sequential adsorption process with hard objects is less and less 
efficient as the dimension increases. (For the RSA of spherical particles, ^(oo) ~ 0.75 for 
D = 1, 9{oo) ~ 0.55 for D = 2, 9{oo) ~ 0.38 for D = 3, etc; it is, however, interesting 
to note that, as a rule of thumb, the saturation coverage in D dimensions is rather well 
approximated by that in one dimension raised to power D 



An additional motivation for studying one-dimensional models is that they are often 
amenable to exact analytical solutions, as illustrated by the case of the car parking problem 
(see section [11 A| ). The underlying reason for solvability is a shielding property 0, that is 
generally found only in one dimension. The shielding property in one dimension means that 
an empty interval of a minimum, finite length separates the substrate into two disconnected 
regions such that adsorption in one region is not affected by the configuration of particles 
in the other region ||^. As a result, a closed system involving only a finite number of kinetic 
equations can be written and, often, solved. Further examples will be given in sections 
and|VI|. 

D. Monte Carlo algorithms 

Computer simulation provides additional information on sequential adsorption processes. 
It was from simulation results, for instance, that Feder first observed the characteristic 



power-law kinetics at long times |^. The basic procedure for a Monte Carlo simulation 
of a sequential adsorption process with hard-core particles is the following: at each step, 
one attempts to place a particle in the simulation cell with periodic boundary conditions, 
at positions sampled from a uniform probability distribution. Because of the finite range, 
hard-core interaction between particles, it is advantageous to employ a grid construction 
T6| . The computer time then depends on the number of particles instead of the square of 
the number. Moreover, close to the jamming limit the process is dominated by very few 
successful events of adsorption and accelerating procedures have been introduced. Indeed, 
at long times, the available surface for inserting new particles is very small compared to the 
simulation cell. If the grid cell element is chosen for accommodating only one particle, the 
algorithm can be improved as follows: first, an unoccupied grid element is chosen at random 
and a second random number determines the trial position within the vacant grid element; 
the time is then incremented by an amount At = N total /N vacant where Ntotai is the total 
number of grid elements in the simulation cell and N^acant is number of vacant grid elements 
at a given density P7| , PB[ . (The procedure accelerates significantly the simulation, but it 



may occur that in the last stage of the process, the surface which is the sum of the grid 
elements remaining vacant exceeds largely the available surface for inserting new particles 
and many unsuccessful attempts still slow down the simulation). Event-driven algorithms 
have been used in various adsorption models (RSA of nonspherical particles [^, polymer 



adsorption |3^). Moreover, in some cases, e.g., in one dimension, the available surface can 
be determined exactly, and a trial position can be always chosen within the available surface 
such that no rejection (or wasted time) occurs along the simulation. The procedure leads 
to a correct result provided that the time is incremented by the ratio of the system size to 
the available surface. It can also be used in the simulation of processes involving desorption 
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(see Section [VII Q ) . 



IV. ADSORPTION OF NONSPHERICAL PARTICLES AND MIXTURES 

A. RSA of anisotropic particles 

Biological molecules have in general a nonspherical shape and they adsorb in an irre- 
versible way when the area of the surface they have in contact with the substrate is maxi- 
mum. Experimentally, Schaaf et. al |^ found that the maximum coverage of the substrate 
they were able to obtain in the adsorption of fibrinogen (a nonspherical protein with aspect 
ratio of approximately 7.5) was around 40%, i.e. less than the saturation coverage predicted 
by the RSA of hard disks (~ 55%) and observed in experiments involving fairly spherical 
globular proteins [||,|1^,0]. (A similar effect was observed for adsorption of albumin [0.) A 



simple way to account for these effects consists in considering the irreversible adsorption of 
anisotropic particles onto a surface with the RSA rules generalized so that the positions and 
the orientations of the adsorbing particles are chosen randomly from uniform distributions 
(for generalizations, see Section |V|. We and other authors p7| , p8| , ^2| -^ have performed this 



study for hard convex bodies of different shapes: ellipses, squares, disco-rectangles (i.e. "2£) 
spherocylinders"). The relevant questions include the following: 

1. What is the influence of the aspect ratio on the saturation coverage of the substrate? 

2. What is the dependence on the particle shape of the kinetics (short and long times)? 

3. What are the similarities and differences between the RSA generated configurations 
and equilibrium configurations at the same surface coverage? 

The hard convex bodies mentioned above are characterized by one parameter which is the 
aspect ratio, a: for ellipses and disco-rectangles, a is the ratio of the long axis to the short 
axis and for rectangles is the ratio of the length to the width. Both disco-rectangles and 
ellipses reduce to a disk when a — > 1, while rectangles tend to a square (i.e. an anisotropic 
object) in the same limit. Fig. |^ displays the saturation coverage obtained by computer 
simulation of two-dimensional RSA as a function of the inverse aspect ratio, 1/a, for various 
shapes. Three main conclusions can be drawn from the results: (i) for the three different 
shapes, the maximum coverage is obtained for a value of a close to 2; (ii) the maximum 
coverage for ellipses and disco-rectangles (for a ^ 2) is larger than that for disks; and (iii), 
for all objects considered here, the saturation coverage goes to zero when the aspect ratio 
a goes to infinity, seemingly with a power law dependence. It is also interesting to note 
that for a ~ 7-8 (typical of fibrinogen, see above), the saturation coverage falls below the 
hard-disk value and is close to 40%. 

As far as the kinetics are concerned, computer simulation also showed that the approach 



toward saturation is slower for anisotropic objects than it is for disks ||2^. From a geometric 



analysis of the available surface, Talbot et al. [27| showed that the long-time kinetics results 



from adsorption of particles in two types of targets, the larger ones being independent of 
the orientation of the trial particle (nonselective targets) and the smaller ones corresponding 
to the insertion of a particle with a small range of possible orientations (selective targets). 
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To minimize size effects while still retaining the effect of anisotropy, Viot and Tarjus ||2^ 



performed extensive simulations of the RSA of unoriented squares, and the results clearly 
showed that the asymptotic approach follows a power law with an exponent equal to 1/3. For 
ellipses and disco-rectangles with moderate and large aspect ratios, the long-time kinetics 
are described by the same power law, but the crossover time beyond which the asymptotic 
kinetics sets in increases with a. For weakly elongated objects, by introducing a small 
anisotropy parameter, e = a — 1, it can be shown that the long-time kinetics follows the 
equation |^ 

p(oo) - p{t) - e (^^^^ + Ce-^M (35) 



where z = e^t and C is a positive constant. Eq. ( ^5|) can be simplified in the two limits 
z << 1 and z » 1. The long-time kinetics thus consists of two consecutive asymptotic 
regimes. First, for 1/y/e >> t >> 1, one has a t"^/^ behavior and for t >> 1/y/e this regime 
vanishes and a t~^/^ behavior emerges. This last regime disappears and the contribution 
to the coverage goes to zero when the aspect ratio goes to one. Simulations have been 
performed for weakly elongated objects and only an effective exponent between 1/2 and 1/3 
has been measured |^. (To observe numerically a clear separation of these two critical 
regimes, simulations of very slightly elongated particles should be performed.) 

The low-coverage expansion of the available surface function $ permits, as shown in 
section [111 A|, a systematic description for low and, somewhat intermediate coverages. To 



third order, one has 



$ = 1 - ^2^ + (2Bl - ^53) e^ + [3B2B3 - ^{Bl + 54) - B'^]e^ (36) 



where Bi denotes the ith. equilibrium virial coefficient I^Of and B'^ is the specific RSA co- 
efficient which can be expressed in terms of Mayer diagrams |^. All coefficients must 
be calculated numerically in general |Q, except the second virial coefficient for which an 
analytical formula is known |^9|. The comparison with simulation results shows that the 



third-order expansion provides a good description as long as $ is not too small. It is worth 
noting that the range of usefulness of the density expansion diminishes when the aspect 
ratio a increases. Even at low coverage, $ decreases markedly when a increases, which 
means that adsorption of very elongated particles rapidly becomes difficult. An interme- 
diate regime appears before the asymptotic regime, in which the longest linear size of the 
particles is responsible for many rejections of trial particles in the process. As displayed in 
Fig. |l|, this regime leads to the formation of a local order whose (linear) domain size is given 
by the length of the particle long axis. The feature becomes more pronounced when the 
aspect ratio becomes larger. This regime cannot be described by a low-density expansion 
whose validity (or quality) range is of order 1/a. (Notice that for equilibrium systems, a 
virial expansion is also unable to describe the isotropic- nematic transition.) To understand 
the emergence of this third regime, it is useful to consider the adsorption of infinitely thin 
needles whose length corresponds to that of the long axis of elongated particles. For needles, 
which do not have a proper area, no saturation is reached and the number density of needles 
increases algebraically with time [HSl . By studying the one-dimensional problem, two of us 
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45| showed that the RSA of needles can be mapped at long times onto a solvable two vari- 



able fragmentation model fl^JSQ] and that the density of needles increases asymptotically 
as 

p{t) ~ t^-\ (37) 

The same algebraic law is expected to hold in two dimension, as is confirmed by computer 



simulations [51, 



Following the procedure described in section [1IIB| , approximation schemes can be con- 
structed by combining low-density expansions and asymptotic regimes. For weakly and 
moderately elongated particles ( a < 5), a fit of the form of Eq. (|3^ gives results which 



are indistinguishable from the simulation data during most of the process. However, the 



predicted saturation coverage is overestimated by the fitting formula [ ^ . 

For very elongated particles, a more efficient approximation scheme can be built. From 
the discussion above, it can be assumed that the filling process is mainly composed of two 
critical regimes: a needle-like regime followed by the true asymptotic regime. Matching the 
two behaviors leads to the following result, 

e^{a) ^ a~'/^'^'^'\ (38) 

and the crossover time between the two regimes is 

t, ~ a3/(i+2v^) (39) 

when a -^ oo. These two results are well supported by simulations: for elongated rectangles. 
Vigil and Ziff [^p2| obtained an algebraic law ^oo(«) ~ l/a°'^''~°'^^, whereas the value of 



the exact exponent of Eq. (p8|) is —0.2612 . . .. The increase in crossover time has also been 



reported but not quantitatively estimated |43 



A structural analysis of the configurations of anisotropic particles generated by RSA has 
been performed by considering the pair correlation function |^. For anisotropic hard-core 



particles in two dimensions, the pair correlation depends on three variables: the distance be- 
tween the two particles r and the orientation 6i of each particle with respect the interparticle 
vector. The center-to-center pair correlation function 5'o,o('"); which is the average over the 
orientations of two particles, i.e. has no orientational dependence, already provides useful 
information on the structure. In particular, the highest peak in 5'o,o('") is for a distance r 
close to 1/2(1 + a) when a < 5, whereas the peak shifts towards r = 1 with a secondary weak 
peak emerging for distances close to a when a > 5 (all distances are measured in units of the 
short axis). These results show that the most probable relative orientation of two particles 
evolves from perpendicular to parallel when the anisotropy parameter becomes larger. Al- 
though the center-to-center pair correlation function is always finite at contact, one observes 
a logarithmic divergence of the surface-to-surface pair correlation function at saturation |^ . 



This is reminiscent of the behavior observed in the RSA of spherical particles. As already 
stressed, nematic order is never observed, regardless of the aspect ratio. 

B. RSA of mixtures 

For a general mixture of n components of arbitrarily hard particles adsorbing irreversibly 
on a planar surface, one would like to be able to predict the time dependent and saturation 
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coverages. Even in one-dimension, this is a challenging problem. However, one can deduce 
the behavior in some limiting cases. 

The RSA of a binary mixture of hard disks is characterized by two parameters: k = 
kA/ks, where k^ and ks are the adsorption rates of the components on an empty surface. 



and the diameter ratio A = aA/o'B- Talbot and Schaaf |^3[ analyzed the case of a binary 
mixture of hard disks of greatly differing diameters, A << 1. The time-dependent coverages 
can be obtained from the numerical solution of two coupled first-order differential equations, 

xbOb) (40) 
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where x^ = (1 + \~^Y.,xb = (1 + A)^ and 6*^ = ^a/[1 — xb^b] is the effective coverage 
of the smaller disks. As long as A is sufficiently small, (pBB depends only on 6b and the 
large disks approach their saturation coverage exponentially. The small disks then behave 
essentially as a one-component system in a reduced area, approaching their jamming limit 
coverage according to the usual algebraic power law, t~^/^ [Q. The final combined saturation 
coverage is accurately given by 

9aMoo) = ^b(oo) + (1 - xb9b{oo))9^ (42) 

where 6^ = 0.547. 

Tarjus and Talbot examined the asymptotic approach to the jamming limit of a contin- 



uous polydisperse mixture ||2^ . The adsorption rate of a particle of radius a on an empty 
surface, K{a), is assumed to be continuous for ai < a < a2 and zero otherwise. As in the 
monocomponent system (section |111B|) , the asymptotic kinetics are determined by the filling 
of isolated target areas. If K{ai) is different from zero then it follows that the jamming limit 
is approached as 

Poo - p{t) ~ t-1/3 (43) 

This is consistent with the idea that the exponent is the inverse of the number of degrees of 
freedom of the adsorbing species: two translational, plus one corresponding to the continuous 
distribution of particle diameters (see also section pVA| for anisotropic particles). Note that if 
K{ai) = then the exponent is determined by the order of the first non- vanishing derivative 



of K at ai 29 



Meakin and JuUien [0,0 conducted computer simulation studies of a binary mixture 
of large and small disks (and spheres), as well as continuous mixtures with uniform and 
gaussian distributions. Their results confirmed the theoretical analyses. 

V. BULK TRANSPORT ISSUES 

In the simple RSA model the position of the trial particle is selected from a uniform, 
random distribution. While this simple choice facilitates both simulation and analytic solu- 
tions in come cases, it is not clear that it is always consistent with the transport mechanism 
of the particle from the bulk to the vicinity of the adsorption surface. 

17 



An adsorbing particle is in general subject to Brownian, gravitational and hydro dynamic 
forces, as well as specific interactions with the adsorbing surface and the pre-adsorbed parti- 
cles. Examples of the latter include van der Waals, electrostatic and short range repulsions. 

A substantial body of work has addressed the issue of the effect of the transport mecha- 
nism of the adsorbing particles on the adsorption kinetics and the resulting structure of the 
deposited particle configurations. In particular, two limiting cases of interest are Diffusion 
Random Sequential Adsorption (DRSA) and the Ballistic Deposition (BD). In the former 
the transport mechanism is pure diffusion, while in the latter gravitational forces are dom- 
inant. For spherical particles, one can define a dimensionless gravity number (proportional 
to the Peclet number) 

where Ap is the difference between the particle and solution mass density, g the acceleration 
due to the gravity, a is the particle diameter, ks is the Boltzmann constant, and T is the 
temperature. The limits Nq ~^ oo and Nq -^ correspond to the BD and DRSA models, 
respectively. 



A. Bulk diffusion 



Schaaf et al. analyzed the effect of diffusion on the asymptotic kinetics |^ . They found 
that the presence of bulk diffusion modifies the asymptotic kinetics, 

Poo - p{t) ~ t-"\ (45) 

The saturation coverage is approached more rapidly than in simple RSA due to a funneling 
effect of the pre-adsorbed particles forming the target. Yet, this power law has never been 
observed experimentally. This may result from the presence of hydrodynamic interactions 
between the diffusing particles and the particles on the surface, interactions that are ne- 
glected in the diffusion RSA model and that may modify the kinetics and the structure of 
the adsorbed layer [p7|-|59|. 



Senger et al H^eI] examined the kinetics and saturation coverage of the DRSA process 
by simulating the random walk of a spherical particle of radius i? on a cubic lattice with a 
lattice parameter a. Each trajectory starts from a randomly selected lattice site in a plane 
at height z = 3R. A particle is adsorbed, and remains permanently fixed, once its center 
reaches the plane z = R. If the particle reaches the plane z = 5R it is considered lost to 
the bulk and a new trajectory is initiated. At each Brownian step, the particle moves to 
one of six neighboring sites with equal probability. If the selected displacement results in 
an overlap with a preadsorbed particle, the particle is returned to its initial position and a 
new direction is chosen. 

The principal result of the simulations was that to within the statistical error, the cover- 
age and structure of the adsorbed particle configurations generated by a DRSA process are 
identical to those of the simple RSA model. In order to gain insight into this unexpected 
result, various one-dimensional models were studied. In particular, the generalized parking 
process proposed by two of us |6^ proved useful in understanding the effect of the transport 
mechanism. 



In (1 + 1) dimensions, the kinetic equation describing the generahzed parking process is 
dG{h,t) 



dt 



POO 

-ka{h){h-a)G{h,t) + 2 ka{h')G{h',t)dh' (46) 

Jh+<T 



Where the rate of adsorption per unit length in a gap of size h is denoted by ka{h). In 
simple RSA, ka{h) = kaH{h — a) where H{x) is the Heaviside unit step function. 

In the generalized parking process, the rate of deposition of a particle in a gap depends on 
the width of the gap, but is uniform within the gap. It has been shown that all generalized 
parking processes have the same jamming limit coverage. Here we reproduce the argument 



given by Bafaluy et al. [Q. The problem is to determine the number, N^{h) of particles 
adsorbed on a line segment of length h after an infinite time. Insertion of one particle into 
this gap leads to two new gaps of length h' and h~ h' — a so one has the following recursion 
formula: 

N^{h) = 1 + 2 r~" N^{h')P{h, h')dh\ (47) 

Jo 

where P(/i, h') is the probability that insertion of a disk into the gap of length h produces 
gaps of length h' and h — h' — a. In all generalized parking processes, including simple RSA, 

P{h,h') = {h-1)-^ (48) 

since, by definition, the trial particle arrives randomly and uniformly in the gap. Thus, 
the final state of the system does not depend on k{h) and we conclude that all generalized 
parking processes have the same jamming limit. The kinetics, of course, can vary greatly 
depending on the form of k{h). 

The story is not yet complete, however, since a careful simulation study revealed that 



the saturation coverage of ID DRSA [^ at 0.7529 is slightly larger than that of ID RSA 
(0.7476 . . .). This difference is a result of the non-uniform flux of particles at the surface. 
Eq. (EUp can be modified to allow for a non-uniform distribution of incoming particles: 



dG{h,t) 
di 



roo 

-k{h)G{h,t) + 2 k{h',h)Gih',t)dh' (49) 

Jh+a 



where k{h) is the total rate that gaps of length h are destroyed by addition of a new particle 
and k{h', h) is the probability per unit length per unit time that deposition of a disk in a 
gap of length h' produces gaps of length h and h' — h — a. One can also obtain information 
about the jammed state without solving Eq. ( ^91) directly with an approach similar to that 
described earlier for the generalized (uniform) parking process. It is possible to obtain an 
analytic expression for the flux of particles on a line segment bounded by two preadsorbed 
disks |6^ and the result obtained for the jamming coverage is equal to the simulated value 
within the confidence interval. More recent work has examined situations where both gravity 



and and diffusion play a role [|65 |. 



B. Ballistic transport and generalizations 

The ballistic deposition model describes situations in which the transport is dominated 



by gravitational effects, corresponding to large values of Nq (Eq. (^J))- Adsorbing particles 
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follow linear trajectories from the bulk to the adsorption surface. If the particle arrives at an 
unoccupied region of the surface, it is immediately accepted and remains permanently fixed 
in place. If it should land on top of a pre-adsorbed particle, it is not immediately rejected 
as with simple RSA. Instead, it follows a path of steepest descent over the pre-adsorbed 
particles. If it reaches the surface, it is accepted, while if it remains in an elevated position 
it is rejected. The possible trajectories for 2+lD BD are shown in Fig. ^ 

1.1 + 1- dimensional models 

Like simple RSA, the model is exactly soluble in one-dimension with a gap-density 
approach [^]. This a consequence of the already mentioned shielding property (see sec- 
tion [111 C| ) . Specifically, the kinetic equation is 



dG{h,t) 

d{kat) 



-H{h-(r){h + a)G{h,t) + 2G{h + a,t) + 2 G{h',t)dh', (50) 

Jh+a 



where ka is the rate of adsorption of disks of size a. A gap of length h will be destroyed if 
the center of an incoming disk falls anywhere on an interval h + a centered on the gap of 
length h. Conversely, a gap of length h may be created by the impact of a disk on either 
of the adsorbed disks bounding a gap of length h + a, or by direct deposition in a gap of 
length h' > h + a. To solve the kinetic equations, one sets G{h,t) = e~''''^'^'^"^^F{kaat)/a'^ 
for h> a, which leads to a solvable differential equation for F{t). G{h,t) for h < a can be 
found using Eq. (|50|). By employing Eq. (^), the solution for the time-dependent density is 
then obtained as 

p(i) = ir°*<i-(l + 2.)e^('-«— ):^ (51) 

where F{t) is given by Eq. (^. As expected, the saturation coverage, p(oo)cr = 0.80865.. 
is greater than that of the RSA process. Moreover, the saturation coverage is approached 
exponentially: 

where 7 is the Euler constant. The saturation state is then obtained more rapidly than 



the corresponding RSA processes (see section |IIA|) . Another distinct feature of BD is the 



presence of connected clusters of particles, which result from the rolling mechanism. 



Viot et al. introduced a model that generalizes both RSA and BD |^. A disk that 
arrives on an unoccupied surface is accepted with probability p. Otherwise, if it alights 
on a pre-adsorbed particle, it follows the path of steepest descent. If this disk reaches the 
surface, it is accepted with probability 1 — p, otherwise it is rejected. RSA is recovered when 
a = p/(l — p) = 0, while a = 1 corresponds to the BD model. As a — ;> 00 only deposition via 
rolling is permitted. This limit corresponds to an Eden-type off-lattice ballistic aggregation 
model. 

In (1 + 1) dimensions the kinetic equations describing the process are 
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^Z,l^ = -Hih -a)ih-a + 2aa)G(h, t) + 2aG(h + a,t) + 2 T dh'G{h\ t), (53) 

d[kat) Jh+cr 

and they may be solved analytically. The number density is 

p(t) = - t''' du{l + 2aM)e2-(i-"-'^-"):^ (54) 

a Jo M 

For kacrt » 1/a, the saturated density is approached as 

-2kacrat 

kaCrH 



p(oo)-p(t)~2e2(-^)V^-, (55) 



showing that saturation is approached exponentially for a > while the usual power law is 
recovered for a = 0. For small, but finite values of a, an intermediate critical regime for 1 << 
kaCrt << a is present in which the number density increases like 1/t; the final exponential 
regime has a contribution that vanishes when a — >■ 0. (This absence of discontinuity, when 
the tuning parameter a decreases to zero, is similar to the situation encountered in the RSA 
of anisotropic particles when the anisotropy parameter goes to zero: see section |I V A| . ) 



The rolling mechanism leads to the formation of connected clusters of particles of different 
sizes, the distribution of which is very sensitive to the value of the tuning parameter a. 
Some examples are shown in Fig. |^. The presence of clusters has a signature in the pair 
correlation function which has an infinity of singularities for each integral multiple of a 
Ip^ . The amplitude of these singularities decreases super-exponentially with the distance. 
Nevertheless, no percolation transition is expected in one dimension for finite values of a. 

2. 2 + 1- dimensional models 

The behavior of the model in (2+1) dimensions is qualitatively similar to the (1+1)D 
version [^,^: a saturation coverage of 0.611 (higher than the RSA value) and an expo- 



nential approach to saturation. As discussed in section |III C| , the irreversible nature of the 
process allows us to determine he asymptotic behavior, to leading order, by investigating the 
filling of targets. Characterizing the targets by the the surplus area relative to the minimum 
value s = S — Sm, where Sm is the smallest target defined by three preadsorbed spheres, 
the density n{s,t) of such targets, when s ^ 0, evolves in the asymptotic regime, t > tc, 
according to 

d 4 

— n(s,t) = -aka-{Sm + s)n{s,t), (56) 

at n 

where one neglects the (much less efficient) filling by direct deposition. One finds that, as 
in the (1 + 1)— dimensional model, the coverage is approached exponentially |[70[| : 

-aka<J^2^t 



Using the methodology outlined in section [III A| and [1IIB| , one can also perform density 



[or coverage) expansions. The generalized ballistic deposition model can indeed be shown 
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to be equivalent to a purely two-dimensional cooperative sequential adsorption model of 
hard disks in which an incoming disk interacts with preadsorbed disks through effective 
potentials involving up to 5— body irreducible interactions. For instance, the 2— and 3— body 
normalized adsorption rates (see Eq. (^Ij)) are expressed as 



e(,ri|r2J 
e(ri|{r2,r3}) 



H{ri2 - a) + a5{ri2 - a) 
e(ri|r2)e(ri|r3) 



+ 2*^(^12 - 0')^(^13 



(y 



sm(ai,23J - 2 



(58) 



(59) 



where if (ri2 — cr) is the Heaviside step-function that describes a hard-disk interaction, 5{ri2 — 



a) is a Dirac function that describes a sticky-disk interaction [^] and 0:1,23 is the angle 
between the vectors ri2 and ri3 (see Fig. ^). Using then the Kirkwood-Salsburg-like hierarchy 
described in Eqs. (^11)- (p9D, expanding in density, and calculating analytically the various 
geometric factors, one can obtain the rate of adsorption $(a, 9) to the third order in coverage: 
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(60) 



Eq. (|60D reduces to that of RSA and of BD for a = and a = 1, [^] respectively. In 
the latter case, the the first- and second-order terms of the expansion vanish, expressing the 
inability of one or two pre-adsorbed particles to prevent an incoming particle from reaching 
the surface. 

Note finally that although configurations of particles generated by the simple ballistic 
deposition model (a = 1) do not contain any spanning clusters, above a threshold value of 
approximately a^ = 3.05 there is a percolation transition as the surface coverage increases 
and one can construct a percolation phase diagram of the model |7^ . We also showed that 
the critical exponents (3, 7, u are consistent with those of ordinary 2D lattice percolation. 

The ballistic deposition model and its generalization are useful to describe the adsorption 
(or deposition) of particles that are denser than the solvent and are therefore subject to 
gravitational forces. This is often the case for colloidal particles. Wojtaszczyk et al. ||2^ 
in their experimental study of the deposition of melamine particles on mica surface showed 
that the saturation coverage ^00 is close to 0.61, and that it is reached with a fast kinetics. 
These features are describable by the two-dimensional ballistic deposition model (a = 1), 
and the connection can even be made more precise by comparing the experimental pair 
correlation function in the deposited layer with the predicted pair correlation function of the 
balhstic deposition model [^ (see Fig|^). For particles not as dense as melamine particles, 
saturation coverages intermediate between the RSA (6*00 — 0.55) and BD (6*00 — 0.61) values 
are obtained |]75[ . Recently, Csiics and Ramsden |7B[ have shown that the adsorption of bee 
venom phospholiphase A2 on a planar metal surface can be described by this model with 
a = 0.47. 
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VI. SURFACE EVENTS: CONFORMATIONAL AND ORIENTATIONAL 

CHANGES 

In the standard RSA model, adsorbing particles are assumed to remain indefinitely in 
their initial adsorbed state. This feature is in contrast with many experimental observations 
of non-spherical and/or flexible proteins and other biomolecules whose adsorbed state may 
change over time. For example, the highly elongated protein fibrinogen (aspect ratio of 
7.5) is known to adsorb initially in an end-on orientation and, following some time, convert 



to a more stable side-on orientation |£3]. Other proteins, upon adsorbing to a surface, 
exhibit changes in secondary or tertiary conformation. A number of early and more recent 
experimental studies have investigated these post-adsorption transitions for various proteins. 
In previous articles, we review these findings in some detail |)75|-p(]|, P^ , pT| ] . Here we simply 
summarize the emerging picture of the transition process: 

1. Post-adsorption transitions in conformation and orientation do occur frequently in 
protein adsorption systems. 

2. The transition usually leads to a larger contact region, thus decreasing the probability 
that incoming proteins land on unoccupied surface. 

3. The transition usually leads to stronger surface-binding, thus decreasing the rate of 
desorption. 

4. The transition tends to be disfavored when the surface is crowded, probably due to 
steric blocking by neighboring proteins. 

It is clear that the presence of a post-adsorption transition affects the kinetics of the 
adsorption process and vice-versa. A truly universal adsorption model should incorporate 
the possibility of a post-adsorption transition, should account for the above experimental 
observations, and should predict, in addition to the total adsorbate surface density, the 
fraction of molecules in an altered state. In the remainder of this section, we review recent 
work aimed at developing RSA-like models that account for a post-adsorption transition. 
The basic idea is that particles first adsorb to the surface as in RSA and then may change 
size subject to certain geometric exclusion rules. The change in size represents a change 
in conformation or orientation. We also highlight the phenomena displayed by these new 
models and the link to experimental results. 

A. One-dimensional models 

Recently, three one- dimensional models of irreversible adsorption with subsequent tran- 
sition have been introduced |j78|. In each of these, particles are modeled as line segments 
of initial length aa that deposit randomly and sequentially onto an infinite line at a rate 
ka- Once placed, the particles may spread immediately to a larger size a^. One of these 
(Model I) incorporates a symmetric spread to the larger particle length. Another (Model 
II) accounts for an asymmetric spread due to contacting another particle during the spread. 
The third (Model III) accounts for a tilting of the particle to either side following adsorption. 
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(See Fig. |^). In each case, the transition occurs only if space allows (i.e. the presence of other 
particles may block the transition). If the transition is assumed to occur instantaneously 
following adsorption, these models become exactly solvable by introducing the gap density 
function G{h,t), defined so that G{h,t)dh is the density of empty line segments of length 
between h and h + dh a.t time t. A kinetic equation can be written for the time evolution of 
G{h, t) for each of these models as follows. 

^, '/ = -{h - aa)G{h, t) + 2 / dh'Gih', t) + 2 / dh'G{h', t) Model I 

o[kat) Jh+aa Jh+ap 

(61) 



dG{h,t) 

d{kat) 



r+oo 

-{h- a^)G{h, t) + 2 / dh'G{h', t) + {ap - ao)G{h + a^, t) Model II (62) 

Jh+aa 



dGih t) /■h+2a3-ac /•+00 

— f^ = -{h- a^)G{K t) + 3 dh'Gih', t) + 2 / 

0\Kat) Jh+afj Jh+2ai3~(7c, 

(63) 



rh+2crg — aa f+oo 

ih-(ra,)G{h,t) + 3 dh'G{h',t) + 2 dh'G{h',t) Model III 

Jh+aa J h+2aa~arv 



In each of these, the first term on the right accounts for the destruction of gaps due to 
incoming particles landing on empty line segments and subsequent terms account for the 
creation of gaps formed by a previously adsorbed particle and one landing on the line (the 
factors of 2 account for the left-right symmetry of the line)Q One can also write similar 
equations for gaps of /i < cTq,, but they are not needed to determine the surface densities. 
To arrive at a solution, one takes G{h, t) = F{t) exp[—ka{h — aa)t] and solves for F(t) (note 
that F(t) is different for each model). 

The kinetic equations for the time evolution of the particle densities, pa and p/j, are 
directly calculable from the gap density distribution via 



d{kj) 

dpp 



d{kat) 



/ dh{Mm{af3, h) - a^)G{h, t) (64) 

J a„ 



dh{h-ap)G{h,t) (65) 



These are integrated to obtain an analytical expression for the densities as functions of time. 
In each model, the line fills initially with /5-particles. At short times, one obtains p/3 ~ t 
and Pa ~ t^ (Model I) and pa ~ t^ (Models II and III). The exponents reflect the number 
of particles on the line needed to block the transition (one in Model I and two in Models II 
and III). Line filling later in the process is dominated by a-particles; these approach their 
saturation as t~^ . /3-particles, conversely, approach saturation exponentially. 



^This property is lost for a class of models in which the incoming direction of particles is not 
vertical (shadow models) [^ 
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In all of these models, the saturation values of the a-particle density and total coverage 
{9 = Pa^^a+Ppc^p) increase and those of the /^-particle density and total density decrease with 
increasing ap/aa- Interestingly, the /^-particle coverage increases with af^/aa in Model II and 
decreases in Models I and III (see Fig. |10^). This quantity is important experimentally since 
often the unaltered fraction may be removed by surfactant elution. The observed behavior 
is due to the greater efficiency of the spreading mechanism in Model II. This spreading 
efficiency is also evident in the average particle diameter, 9/ p. (see Fig. |TO|b). 

B. Two-dimensional models 

The experimental situation of interest is usually adsorption onto a two-dimensional sur- 
face. With this in mind, recent work has focused on the development of two-dimensional 
models of irreversible adsorption with a post-adsorption transition [^,|^,^,|^ . Adsorbing 
particle are modeled as disks of size cTq that adsorb randomly and sequentially onto a plane 
at a rate ka = kaC (c is the bulk solution concentration). Once adsorbed, if space is available, 
the disk will expand symmetrically and instantaneously to a larger diameter a^ at a rate kg- 



If space is not available, the disk remains permanently of size aa (see Fig. |Tl]). The kinetic 
equations for the time evolution of a- and /3-particles are 

Op 

-—^ = kaC^aip, S, Ks) - ksPa'^alsiP, S, K^) (66) 

dpp 



dt 



hPa^aM^^^s) (67) 



where $„ is the probability that an a-particle lands on unoccupied surface and "^a/s is the 
probability that sufficient space exists for the transition to occur. These functions depend 
on the overall surface density, p, the ratio of particle diameters, E = (T/s/cTa, and the ratio 
of the spreading to adsorption rates. Kg = kg/kacaa (eta is the area of an a-particle). In 
the special case where ks=0, these equations reduce to those of the standard RSA problem. 
When kg approaches infinity, which is just the two-dimensional analog of Model I above, one 
has kinetic equations 



dpa 
dt 



:M<l>„(p,S)-<l>^(p,S)) (68) 

dpp 



dt 



kac^pip,^) (69) 



where ^p is the probability that an incoming particle will have sufficient space to immediately 
spread to a diameter o"/?. 

Analytical solutions to Eqs. (|6^-p9|) are not available. However, $q,,$^, and \1/q,/j may 



be expressed as power series in terms of the total surface density. This is done by first 
generalizing the Kirkwood-Salsburg-like hierarchy (see Eq. (^) ) to express sets of mixed 
particle-cavity distribution functions in terms of an infinite series of integrals over Mayer 
functions and multiplet density distribution functions (see section [III A|) . 
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*AiA2...A„+„(ri, r2, ..., r„, r^.^^, ..., r;_^_^) is the probability density of finding a Ai-particle 
at position ri, ... a A„-particle at position r„, sufficient space (a cavity) for a A„+i- 
particle at position r„+i, ..., and sufficient space for a A„+m-particle at position Vn+m- 
PA"A2...a('^i'J'2,---,rn)^AiA2...A„_ia/3(ri,r2,...,r„_i,r;) is the probabihty density of finding a 
Ai-particle at position ri, ..., a A„_i-particle at position r„_i, and an a-particle at position 
r„ that has sufficient space to spread and become a /3-particle. The terms in these sum- 
mations represent the blockage of cavities on the surface by Sq, a- and sp /5-particles. All 
of the distribution functions (the $'s, \l/'s, and p's) may be written as series expansions 
in the total singlet density p. Using Eqs. (|70D - (|7TD and kinetic equations, Eqs. (|66D- (|67D 
for the time evolution of the multiplet density distribution functions, the coefficients of the 
density expansion may be determined term by term |8y]. In the case when Kg is finite, a 
renormalization is performed to change variables from (p. Kg) to (p, () where ( = pKg. This 
variable change is reflected in the form of the kinetic equations for all density functions and 
results in the series expansion coefficients having a dependence on the new variable C,. In the 
case where Kg is infinite, no renormalization is needed and the series expansion coefficients 
depend only on S. 

Power series are accurate only at low to moderate densities. As the system approaches 
saturation, one must use other tools to analyze the kinetics. As is the case in standard 
RSA, it is useful to consider the geometry of small, isolated targets on the crowded surface 
that are available for addition of a- and /5-particles. By expressing the rate of adsorption 
in terms of the density and size of the targets, one finds that a-particles approach their 
saturation as t~^/^ and that /5-particles approach their saturation exponentially in time. 
An exception is when S=l; in this case the /3-particles approach their saturation as t~^/^ 
since their formation is not inhibited by surface blockage. Consideration of surface targets 
also allows one to show that the S derivative of the saturation a- and /3-particle (but not 
the total) density diverges as S approaches 1. This result suggests that even a very small 
transition could result in quite a large fraction of particles remaining in the unaltered state. 
Finally, an asymptotic analysis of the structure of the adsorbed layer shows that the (3 — (5 
pair correlation function, gp^i^r), diverges at contact (as does g(r) in standard RSA) while 
gaa{r) and ga/si^) remain finite. 

The asymptotic regime scaling laws can be combined with the short/moderate time 
regime density series expansion via interpolation yielding expressions that are accurate over 
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the entire filling process. When Kg is infinite, methods similar to those developed for 
standard RSA may be applied. One obtains equations for the time evolution of the total 
density and total fractional surface coverage, 

df>_ ^ kgcjl-xf 

dt l + bix + b2X^ ^ ' 

de_ _ kgcapjl -xf 

dt 1 + CiX + C2X'^ + C3X^ 

where x = p/poo and the coefficients 61, 62, Ci, C2, C3, and poo are determined by matching the 
density expansion of the Pade ratios with the coefficients determined above. These equations 
may be integrated to give pa and p^ as functions of time that compare quite favorably with 



computer simulation 80 



When Kg is finite, it is more accurate to introduce an effective area, ae//(p, C), that 
serves as a mapping variable of the RSA+spreading model onto the standard RSA model. 
It is defined as the particle area in the standard RSA of disks that would yield the same 
fractional available surface as in the RSA+spreading problem at the same overall surface 
density, that is, $(6* = paeff) = ^o(p) where $ is the available surface function for the RSA 
of disks. Also introducing the average area of particles landing on the surface of density p, 
a(p, 0, kinetic equations for the density and coverage may be written: 

do 

-^ = kac^pa,ff{p,0) (74) 

— = kaCaf3^pa{pX)) (75) 

Qeff and a are themselves written as Pade approximants whose coefficients are determined by 



matching the low density expansion coefficients |^ . This method yields extremely accurate 
total surface densities and reasonably accurate partial densities when compared to computer 
simulation. 

When Kg is infinite, results of computer simulation of the two-dimensional problem are 
qualitatively similar to those of the one-dimensional problem. One difference is that the 
density of a-particles (/5-particles) tends to be higher (lower) in 2-D due to the greater 
probability of landing near to a preadsorbed particle. When Kg is finite, some important 
differences emerge. First, the a-particle density may become non-monotonic in time due to 
the delayed spreading (see Fig. [1^). Second, the /3-particle coverage, Pjsa^, approaches zero 
for large S. This is because for any finite rate of spreading, one can choose a spreading 
magnitude larger than the characteristic separation of adsorbed particles. Increasing Kg 
favors the spreading event and causes a decrease in the pa and p and an increase in p^ and 

e. 

The need to employ a concentration-dependent particle size when fitting experimental 
data to the standard RSA model |T^ provided an incentive to develop irreversible adsorption 
models incorporating a post-adsorption transition. With the RSA+spreading model, several 
kinetic isotherms, each differing only in bulk concentration, may be accurately predicted with 
a single particle area that is close in value to the known size of the adsorbing particle (see 

Fig. p m. 
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VII. INFLUENCE OF DESORPTION PROCESSES 

Up to this point, we have considered adsorption processes that are strictly irreversible, 
that is, ones in which proteins or other macromolecules remain forever on the surface follow- 
ing adsorption. In many experimental situations, proteins do indeed adsorb tenaciously and 
little or no desorption occurs. However, other examples exist where desorption may occur 
spontaneously or in response to a change in experimental conditions (pH, ionic strength, 
surfactant or detergent concentration). It is clear that, in order to be of general applicabil- 
ity, a protein adsorption model must be able to account for the possibility of desorption. In 
this section, we introduce several models that extend the irreversible adsorption approaches 
presented above to include desorption. 

A. Multistep adsorption / Depletion model 

Protein adsorption need not always be conducted as a single step process. Adsorption 
may be interrupted by a buffer rinse or by exposure to a surfactant or detergent solution. 
This usually results in removal of a fraction of the adsorbed proteins. Due to the non 
equilibrium nature of protein adsorption, we suspect that when adsorption is resumed, the 
behavior may differ nontrivially from an uninterrupted process. 

We have investigated the effect of a desorption (or depletion) step, separating two irre- 
versible adsorption steps, on the structure, kinetics, and saturation density of adsorbed disks 
83| , |8^ . A starting point is to consider the evolution of the density distribution functions 



during a random depletion from an initial density pi as follows: 

— = -fcdnp^"^(ri,r2,...,r„;pi,p) (76) 

It is straightforward to show that this leads to a conservation of distribution functions, 
that is 5f(")(ri,r2,...,r„;pi,p) = ^(("^(ri, rg, ..., r„; pi,pi) for all P2 < p < Pi where pi and 
P2 are the densities before and after depletion ((7*^"^ = p^^^p"-) |8^. This result holds 



for depletion of any immobile collection of particles, including those at equilibrium. An 
interesting implication is that one could create a low-density configuration that would have 
the structure of one taken at a much higher density. This result can be used to analyze 
the kinetics and saturation density of a second adsorption step beginning with the depleted 
configuration. To third order, 

$(pi,P2,p) = $^^^(p) + A32P^(P1 - P2) (77) 

where A32 is the contribution of the third order coefficient of ^^^^ from the pair density. 
8^ Since A32 > |jl2[ , the available surface function for the multistep process is greater than 



that for the simple RSA process and an enhancement in the saturation density is expected. 
Using an interpolation scheme, one can calculate the saturation density and show that it 
is always enhanced compared to the uninterrupted adsorption process and that the maximum 
enhancement occurs when p2 = |pi |Q. An interesting aspect of these results is that in 
order to reach a certain surface density as rapidly as possible, it is sometimes most efficient 
to incorporate a desorption step. This is shown in Fig. H. 
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B. Partial reversibility 

A number of experimental results suggest that in some instances protein adsorption is 
partially reversible, that is, a fraction of the adsorbed molecules may be removed by a buffer 
rinse while others cling to the surface irreversibly |^. These observations are consistent 
with the two-state adsorption picture first put forth by Soderquist and Walton where pro- 
teins adsorb initially in a reversible manner and then later become irreversibly adsorbed by 
changing conformation and/or orientation |^ . We proposed a model of partially reversible 



protein adsorption [^ incorporating Morrisey's observation that neighboring proteins can 
block this transition | |86| | , 

In this model, proteins are modeled as disks of diameter o"q that adsorb onto the surface 
sequentially and without overlap at a rate ka = kaC. Once adsorbed, the proteins desorb at 
rate kd or spread to a larger diameter a 13 at a rate kg (the latter is subject to no overlap 
with other adsorbed proteins). The rate laws are as follows: 

dp ■^ 

-^ = kaC^aiP, ^,Ks, Kd) - ksPa^ap{p, ^, Kg, Kd) - kdPa (78) 

^ = kgp^^i/^^{p,E,Ks,Kd) (79) 

where $„ and "^a/s are, respectively, the probabilities of adsorbing and spreading without 
overlap, T. = ^, Ks = r^^—, and Kd = r^^ where a„ = nal/A. 

^' CTc ' * kaCaa kaCa^ °^' 

Because the adsorption is not strictly reversible in this model, a density expansion ap- 
proach as introduced in sections [II B| and |III B| is not likely to succeed. However, the asymp- 



totic regime may be analyzed by again considering the time evolution of isolated targets on 
the surface for a- and /3-particles. Using this approach, one can show that both pa and pp 
approach saturation as t~^/^ [^|. The change from a faster exponential approach to a slower 



algebraic approach for the /5-particles when Kd becomes non-zero is because, in this case, 
a-particles may spread at long times when cavities form due to desorption of neighboring 
a-particles. (In the special case of S = 1, no steric hindrance to spreading is present and 
an even faster t~^/^ approach is found.) 

To evaluate the above kinetic equations over all times, computer simulation is required 



37| . An interesting result is that for certain parameter values, both the total density and 
the a-particle density may become non-monotonic in time. This is due to the gradual 
replacement of initially adsorbed a-particles by a smaller number of larger /3-particles. This 
behavior is similar to what is observed in experimental "overshoots" where the total density 
decreases at long times |^. Another interesting result is that saturation densities tend to 



be higher than in purely irreversible adsorption. This is due to a more liquid-like particle 
structure, as evidenced by the secondary peaks in the radial distribution function, which 
leads to more efficient packing on the surface. 

The chief utility of this model is its applicability to experimental situations in which 
desorption and post-adsorption transitions occur. Good quantitative predictions of the 
adsorption kinetics and the reversible fraction of lysozyme on silica are obtained (see Fig. |15[) 

MM- 
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C. Parking lot model 

In the adsorption-desorption model, particles are placed in a Z)-diniensional space ac- 
cording to the RSA rules at a constant rate ka and in which an additive uniform desorption 
process takes place such that all objects in the system are subject to removal (desorption) 
at random with a constant rate k^. In the parking lot model, the substrate is a line and 
the objects are hard rods. This one- dimensional model has been solved in some limiting 
cases. When kd = 0, the adsorption is totally irreversible and the process corresponds to 
the Car Parking Problem. When ka = 0, starting with any configuration of particles, the 
process corresponds to the desorption model studied in section |V11A| . The limit k^ — > 0+ 
allows a slow but eventual rearrangement of particles on the line leading to a final coverage 
of unity. (It is worth noting the finite discontinuity of the final saturation density between 
the case kd — »■ 0+ and kd = 0, i.e., RSA). Moreover, the final density is independent of 
the initial configuration of particles on the line; it depends only on the ratio of desorption 
and adsorption rates. For kd -^ 0, accurate descriptions have been obtained p8|j89|] . In this 
case, the process cleanly divides into two sub-processes. The initial phase consists of an ir- 
reversible adsorption and is followed by an infinite sequence of desorption-adsorption events 
in which a rod detaches from the surface and the gap that is created is immediately filled 
by one or two new rods. The latter possibility causes the system to evolve continuously to 
the close-packed state with pa = 1, as @,^] 1 — crp(t) ^ l/ln(t) where t now represents a 
rescaled time. For the general case, where both ka and kd are non zero, a complete solution 
is not yet available. 

The properties of the parking lot model depend only on the ratio K = ka/kd- Large 
values of K correspond to small desorption rates. With an appropriate rescaled time, the 
kinetics is given by 

'-^^m-^i, (SO) 

where $(t), the insertion probability at time t (or density p), is the fraction of the substrate 
that is available for the insertion of a new particle. The presence of a relaxation mechanism, 
even infinitesimally small, implies that the system eventually reaches a steady state that 
corresponds to an equilibrium configuration of hard particles with pcqO" = -ft'$cq(PcqO"), where 
Poq denotes the equilibrium density. At equilibrium, the insertion probability is given exactly 
by 

<l>eq(p) = (1 - pa) exp(-p(T/(l - pa)). (81) 

Inserting Eq. ( PD in Eq. (^) leads to the following expression for the equilibrium density: 

where Ly/[x) (the Lambert- W function) is the solution of x = ye^ . In the limit of small 
K., the isotherm takes the Langmuir form i^PeqO ~ -^/(l + ^)) while for large ii", peqa ~ 
1 — l/ln(i^). At small values of K^ equilibrium is rapidly obtained, but at large values 
the densification displays a dramatic slowing down. An adiabatic (mean-field) treatment 
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consists of assuming that, at any density p(t), the structure of the adsorbate follows an 
equilibrium form. This means that $ is similar to Eq. ( ^T]) with p{t) in place of peg- Denoting 
Sp(t) = pit) — Poo, with Poo = Peq{K), ouc performs a first order density expansion of Eq. ( |80[ ) 



^5p = -rMF{K)6p + 0{6p') (83) 

with 

Tmf{K) ~ IniKf/K when K is large (84) 

equivalent to a relaxation time of K/ In(i^)^ for large K pO]. However, a careful analysis of 



the kinetic equations for the gap distribution function has been done ||9^ and it has shown 
that the relaxation rate F is actually given by 

r-2^^ + o(^-^) (85) 

and thus much smaller than the mean-field prediction, Eq. (0); simulation results are 
compared with both the mean field result, and with the prediction from the gap distribution 
approach, Eq. (^), in Fig. |16[ (It is worth noting that the difference between a mean-field 
treatment |8^ and a more detailed description leads to a difference which only occurs in 
subdominant terms of the kinetics). In addition to the anomalous relaxation behavior of this 



model, the power spectrum of fluctuations exhibits two frequency scales [91] which could 



be a general feature of models with very slow relaxation. Similar power spectra have been 



observed in densification of monodisperse glass beads during a long series of vibrations [93 . 



D. Scaled Particle Theory 

The Parking Lot Model described in the previous section shows both that when the 
desorption rate is very small, the relaxation kinetics is very slow and the whole process 
is not well described by a mean-field approximation. On the other hand, if desorption is 
not too small, an adiabatic approximation provides a fair description of the process. The 
final state is an equilibrium one and the difference between the transient states and the 
corresponding equilibrium states at the same coverage is small. This property is useful 
when one considers the adsorption of complex systems, such as mixtures, because it allows 
a drastic simplification of the analysis by using an approximation which essentially requires 
the knowledge of equilibrium properties. The Scaled Particle Theory (SPT) introduced by 



Reiss et al. [£J] is particularly convenient since it yields simple analytic expressions for the 
available surface functions for systems at equilibrium. 
The kinetic equation generalized to mixtures is 

-^ = ka,icM{e,}) - kdA (86) 

where $j({^j}) is the available surface function of species i. To apply SPT one first approx- 
imates $ by $£9 and then the connection to the equation of state of the equilibrium mixture 



is provided by the relation [21 
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ln$f 



eg 



-^^f/kT 



(87) 



where /if is the residual chemical potential of component i. 

For a two-component mixture of hard disks of diameter cxi and (T2 the explicit equations 
are: 



$f = (1 - 01 - ^2) exp 
where 7 = cr2/cri is the size ratio. 



$^'' = (l-0i-02)exp 



3^1+7-^7-^ + 2)^2 

1-^1-^2 



3^2 + 7(7 + 2)^1 



(1 - 01 - 02)2 



?2 + 7^l)' 



(88) 



1-01-02 



;i - 01 - 02 



(89) 



The time-dependent coverages obtained by integrating Eq. (|^) with Eqs. (^) and (|89D 
agree very well with numerical simulations of the adsorption-process. For certain parameter 
values, the model exhibits an "overshoot" in the time dependent coverage - a phenomenon 
that has been observed experimentally. 

Adsorbed configurations of a single species that can adsorb in two or more orientations 
with respect to the surface can be considered as a two-dimensional mixture. Kinetic and 
isotherms based on SPT have been developed to describe this situation P5|-P^]. The equa- 
tions presented above can be readily generalized to continuous multicomponent mixtures 
and some interesting kinetic phenomena have been reported |S^. Finally, kinetic equations 
based on SPT have also been developed for the spreading model [p9|. 



VIII. MULTILAYER FORMATION 

The rules of the RSA model and its extensions discussed thus far allow only for a single 
adsorbed layer to be placed on a surface. In fact, a number of experimental situations 
exist where adsorption does not stop at one layer but instead continues with later arriving 
particles adhering to previously adsorbed ones. The fabrication of thin film materials via 
solution deposition of macromolecular building blocks is an important example of multilayer 
adsorption and is currently used to make sensors, superconducting materials, and cathode 
ray tubes. The physical, electrical, and optical properties of these thin film materials depends 
strongly on their structure and density. Models are thus needed to predict thin film structure 
and density in terms of the particle properties and deposition conditions. In this section, 
we discuss multilayer adsorption models developed for this purpose. 

A number of lattice and continuum studies of multilayer deposition have appeared re- 
cently. [[7|, |100| , piIT| . Our work focuses on the first layers (typically less than 20 — 30) and 
that it incorporates all three events that can occur when a particle contacts a growing in- 
terface: it may stick to the interface (i.e. adsorb), it may diffuse away from the interface 
(i.e. desorb), or it may remain in contact while descending further toward the surface (i.e. 
roll) [|102|J103|| . 

Particles are modeled as {D + l)-dimensional spheres that descend vertically onto a 
D— dimensional substrate. If a particle lands on the surface, it remains there permanently. 
If it lands on a previously placed particle, an overhang amount 5 is calculated either with 
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respect to the surface or to the contacting particle (see Fig. [VII]] ). The overhang extent 



is compared to two parameters, an adsorption parameter ^i and a rolling parameter 82- If 
(5 < ^1, then the particle will adsorb and remain permanently in its place. \i 5i < 5 < ^2, 
then the particle will desorb. If 5 > 82-, then the particle will roll toward the surface and, if 
its path is unblocked, will adsorb to the surface. If the path is blocked by another particle, 
then the rolling particle desorbs. The rationale behind these deposition rules is that particles 
landing directly on top of other particles will tend to adhere more strongly due to a greater 
contact area. By adjusting the values of 5i and ^2, one can perhaps mimic most experimental 
situations. 

In one dimension, an analytical solution is available for the lowest (surface contacting) 
layer when overhangs are calculated with respect to the surface. This is of particular im- 
portance since in certain experimental situations, higher layer particles can be removed by 
rinsing and changing the solution conditions. To arrive at the solution, we introduce the 
gap distribution function, G{h,t), defined so that G{h,t)dh is the density of gaps of length 
between h and h + dh at time t. The time evolution of this function obeys the following: 

dG(h t) /■^+'5i 

J, \, =~{h + a + 2(^1 - S2))G{h, t) + 2 / G{h', t)dh' 

0[kat) Jh 

+2 / Gih', t) + 2ka{(J - 52)G{h + (7, t) (90) 

Jh+u 

where ka is the rate of deposition and a is the particle diameter. By taking G{h, t) = 
^-kaih+a+2{Si~S2))t pi^j^^^-l-y ^2 ^ -j-j^g integrodifferential equation reduces to an ordinary differ- 
ential equation in t. Its solution is then used to find the surface density via 



dp 



d{kat) 



/ {h + (r-26i)G{h,t)dh (91) 



For higher layer particles and for particle overhang rules, computer simulation is required. 
The simulation procedure follows the above deposition rules and involves a line of length 
256 cr (1-D) or a plane of size 256a x 256a (2-D). 

We find that the saturation density of the lowest layer increases with decreasing 61 or 
62- In the former case, screening of the surface by higher layer particles is reduced. In the 
latter case, the rolling mechanism (a very efficient way to fill the surface) is enhanced due 
to the screening effect. An interesting result is that the radial distribution function of the 
lowest layer shows only a weak correlation between particles and decays to zero a.t r = 2a 
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When either Si or S2 is unity, the use of surface or particle overhang rules leads to the same 
result ||104||. Otherwise, particle overhang rules always result in a more dense lowest layer. 



1-D and 2-D results are qualitatively similar. The density profile of higher layer particles 
exhibits density oscillations whose form may vary markedly depending on the values of the 
parameters used. For example, when 61 = 62 = 1, the oscillations decay within 4 particle 
diameters of the surface [ p.04| | whereas for 61 = .75, 62 = 1, the oscillations continue past a 
height of more than 8 particle diameters (Fig. |T8|). 
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TABLES 

TABLE L Saturation coverage for RSA processes. Estimates from approximants and simula- 
tion results. 



Estimate 



Computer Simulation 



D = 2 hard disks 

D = 3 hard spheres 

D = 2 parallel hard squares 

D = 3 parallel hard cubes 



0.548-0.553 



0.365 [105] 



0.547 [y] 
0.382 1T0^J109| 



0.5625 |3|,|4| 



0.5620 [|07| 



0.40-0.46 H 



0.42-0.43 |108|J109|] 
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FIGURES 

FIG. 1. Typical RSA (a) and equilibrium (b) configurations of disco-rectangles with an aspect 
ratio a = 15 at the RSA saturation coverage 9 = 0.45. Note that the orientational order is purely 
local for RSA whereas a nematic phase is present at equilibrium. 






(b) 



(c) 



FIG. 2. (a) Typical RSA configurations of hard disks for three values of the coverage (a) 
6 = 0.1, (b) 9 = 0.3, (c) 9 = 0.5. The black disks are the adsorbed hard disks and the hatched region 
represents each disk's exclusion area. The remaining white region corresponds to the surface that 
is available for adsorption of a new disk. In (c), close to the jamming limit, the remaining available 
surface consists only of a small "target" indicated by the arrow. 
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FIG. 3. Typical "target" in the asymptotic regime of the RSA of disks close to saturation, (see 
Fig. ^c). Only a small surface area characterized by a linear length h is available for adsorption of 
a new disk (such a disk is shown by the dotted circle). 



D = 2 



= 0.3 




jWvM'A^A''M'-'%"'-'-'-'^^A 




D = 2 



= 0.5 



FIG. 4. Pair correlation function g{r) for the RSA configurations of hard disks for two values of 
the coverage: (a) = 0.3 (b) = 0.5. The continuous line is computed from the Percus-Yevick-like 
integral equation for RSA, and the wavy line represents the simulation result. 
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FIG. 5. RSA of rectangles, ellipses and disco-rectangles onto a plane: saturation coverage, 
9a{oo), as a function of the inverse of the aspect ratio a (for 1 < a < 15). The lines are drawn for 
visual guidance. 




FIG. 6. Generalized ballistic deposition: schematic illustration of the four possibilities for 
inserting particle 3 in the vicinity of two previously adsorbed particles 1 and 2: (a) corresponds to 
a direct deposition, (b) and (c) correspond to a deposition after rolling on one preadsorbed sphere, 
and (d) to a deposition after rolling on two spheres. This process can be described by a CSA of 
disks with effective rates of adsorption (see Eq. (|59D). 
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FIG. 7. Density of clusters formed by s particles at contact ps{t) as a function of time for 
different values of a > and s. (a) a = 0.1, s = 1, 2,3 (b) a = 1, s = 1,2,3,4,5, (c) a = 3, 
s = 1, 2, 3, 4, 5. Note that the density of monomers pi{t) has a maximum at a finite value of time 
for a = 1 and 3 but not for a = 0.1. For o = 3, the efficiency of the rolling mechanism leads to a 
density of dimers larger than the density of monomers at the jamming limit. 
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FIG. 8. Pair correlation function for deposited colloid layer at four different coverages: (a) 
9 = 0.075, (b) 6 = 0.30, (c) 6 = 0.46, (d) 6 = 0.53. Diamonds correspond to the experimental 
results (melamine latex spheres on mica) and the lines are determined from simulation of the 
ballistic deposition model. 
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FIG. 9. Illustration of a two-state transition on the substrate. The particle is a hard rod of 
length Ga when it reaches the surface, (a) Model I: if the neighbors are sufficiently far, the particle 
grows symmetrically up to a length a^, (b) Model II: growth can also occur close to a neighboring 
particle, but the spreading is asymmetric (c) Model III: the particle adsorbs side on and tilts to 
the left or the right if space is available. 
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FIG. 10. (a) Dependence on ap/aa of the saturation coverage of the largest particles, 
6/3(00) = p/s{oo)a/3, for the three models, (b) Time evolution of the mean size of the adsorbed 
particles for the three models and for ap/aa = 1.5, 2.0. 



a-particle 



p-particle 




41 



FIG. 11. Protein adsorption with surface-induced configurational change. Illustration of the 
two-dimensional symmetric spreading model. Once adsorbed, an a— particle spreads symmetrically 
at a rate ks on the substrate if neighboring preadsorbed particles do not prevent the transition. 
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FIG. 12. Kinetics of RSA+spreading model. Effective area interpolation calculation (dashed 
lines) of the reduced total and partial densities as function of reduced time for S = 1.05 and (a) 
Kg = 1 and (b) Ks = 10. Also shown are simulation results (solid lines). 
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FIG. 13. Surface density of adsorbed fibronectin F as a function of time. Curves corresponding 
to bulk protein concentration of 0.05, 0.01, 0.005, 0.002 g/1 are shown. Also displayed are the 
curves predicted from Eq. ( [7^ with ka = I x W^^cm/s and aa = 20nm. More details are given in 
Ref. |I|. 
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FIG. 14. Influence of a depletion step on the adsorption kinetics. The dashed line corresponds 
to a single-step RSA and the full line to two RSA steps separated by a depletion step (with kd = ka). 
In both cases, the process is stopped as soon as a coverage of 0.5 is reached. 
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FIG. 15. A comparison between the model of partially reversible adsorption and the lysozyme 
adsorption experiments of Walgren et al [^. Simulated total (solid line) and /3— particle (dashed 
line) densities are shown along with the experimental total amount for two bulk protein concen- 
trations. (For the determination of the parameters, see Ref. [^). We note good agreement in 
the total amount adsorbed before rinsing (at 3600s). The amount remaining after rinsing is well 
approximated by the /?— particle density at the time (indicated by squares). 
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FIG. 16. Relaxation rate F for the final exponential approach to equilibrium in the Parking 
Lot model. Upper curve: mean- field approximation, Eq. (B^). Intermediate curve: Eq. (p5|) Open 
circles: simulation results. 
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o § (surface overhang) 
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81 < 8 => Adsorption 81 < 8 < 82 => Desorption 8>82 => Rolling 




FIG. 17. (a) Depiction of the surface and particle overhangs, (b) Schematic of the three possible 
events that can occur when an incoming particle contacts a previously placed particle. 



47 



Q. 




Cl 



1.0 

0.8 
0.6 
0.4 
0.2 
0.0 



0.0 




2.0 



6=1.0 
5,=1.0 



'*^''t'''''**''>'-'Mf/t*'H^k^^ 



4.0 
z/o 



6.0 



8.0 



FIG. 18. Density of higher-layer particles as a function of height above the surface in two 
dimensions for different values of 6i and 52- 
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